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ABSTRACT 

Knowledge of when critical crop stages occur and how the environ- 
ment affects them should provide useful information for crop management 
decisions and crop production models. This research evaluated two 
sources of data for predicting dates of silking and physiological matur- 
ity of corn (Zea mays L.). Initial evaluations were conducted using 
data of an adapted corn hybrid grown on a Typic Agriaquoll at the Purdue 
University Agronomy Farm from 1979 to 1981. The second phase extended 
the analyses to large areas using data acquired by the Statistical 
Reporting Service of USDA for crop reporting districts (CRD) in Indiana 
and Iowa from 1969 to 1980. Several thermal models were compared to 
calendar days for predicting dates of silking and physiological matur- 
ity. Mixed models which used a combination of thermal units to predict 
silking and days after silking to predict physiological maturity were 
also evaluated. At the Agronomy Farm the models were calibrated and 
tested on the same data. For each CRD the models were calibrated using 
4 or 5 years of data and tested using 7 different years of data. 

The thermal models were significantly less biased and more accurate 
than calendar days for predicting dates of silking. Differences among 
the thermal models were small. Significant improvements in both bias 
and accuracy were observed when the mixed models were used to predict 
dates of . physiological maturity. The results indicate that statistical 
data for CRD can be used to evaluate models developed at agricultural 
experiment stations. 


INTRODUCTION 

Crop development, or ontogeny, involves complex physiological and 
biochemical processes which are influenced by the crop's environment in 
ways that are still inadequately understood. Temperature and photoper- 
iod are the principal environmental variables which influence develop- 
ment of crops. In some situations, the availability of moisture and 
nutrients also may affect crop development. 

During the past century numerous models to describe the ontogeny of 
various crops as a function of environmental variables, particularly 
temperature, have been proposed. There are many different methods of 
calculating and accumulating temperature or thermal units for corn ( Zea 
mays L.); for example, Cross and Zuber (7) reported on 22 methods for 
corn. The simplest and most broadly researched method is Growing Degree 
Units (GDU). A base temperature for growth of 10°C is subtracted from 
the mean air temperature to give the daily GDU. Modifications of this 
simple method frequently impose some upper and lower limits on the daily 
temperature inputs (4,7,12,19), while other methods consider day and 
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night temperatures separately (5). For corn these limits commonly are 
30° C for the maximum temperature and 10° C for the minimum temperature. 
A GDU index is obtained by summing the daily GDU from planting to the 
stage of crop development desired, usually silking or physiological 
maturity. 

Considerable effort has been directed at trying to predict flower- 
ing and physiological maturity dates of various crops on the basis of 
temperature data. When cumulative thermal units were used to compare 
maturation of corn hybrids at different locations, those with a base of 
1 0° C more effectively described crop development than calendar days (2). 
Gilmore and Rogers (8) studied the development of 10 hybrids and 10 
inbred lines of corn using 15 different methods of calculating thermal 
units. Thermal units calculated using temperatures taken at 3-hour 
intervals did not estimate silking significantly better than those cal- 
culated using daily maximum and minimum temperatures. Differences among 
hybrids in the rate of development based on accumulated thermal units to 
silking were noted. Other researchers also have observed differences in 
rate of development among hybrids (14,15). 

Numerous empirical and theoretical methods of estimating the silk- 
ing and physiological maturity stages of corn have been devised and com- 
pared (1,3,7,10,12,19). Although differences among the methods for 
estimating a particular stage of development were generally small, all 
methods of accumulating thermal units were better indicators of crop 
development than calendar days. 

Stages of development can be estimated very well for corn hybrids 
of different maturity classes using the simple GDU system with a base 
temperature of 10° C (13). Frequent and detailed data on stages of 
development result in better measures of the relationship between crop 
development and GDU than has been indicated by previous studies using 
only one or two stages of development (13). 

The thermal unit accumulation concept assumes that photo period does 
not influence the rate of crop development (19). For domesticated crops 
grown in areas where they are adapted, development may seem to be inde- 
pendent of photoperiod. This is because the photoperiod is either lon- 
ger or shorter than the optimum photoperiod or because the crop is rela- 
tively insensitive to photoperiod. Corn development is influenced by 
photoperiod (1,6). Decreasing photoperiods hasten flowering (i.e., 
silking) and reduce the number of leaves per plant in corn (1). 
Increasing temperatures also hasten flowering but increase the number of 
leaves per plant (1). For corn grown in U.S. Corn Belt, the changes in 
photoperiod are confounded with changes in temperature and are nearly 
impossible to separate in field experiments. Coligado and Brown (6) 
developed a model incorporating temperature, photoperiod, and genetic 
factors to predict tassel initiation of corn. Their model appears to be 
sound theoretically but needs further research to extend it to other 
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stages of development. Although temperatures and photoperiod interact 
to influence the development of corn, particularly tassel and ear initi- 
ation, thermal models are generally accepted as adequate to predict 
growth and development of corn (11). 

In summary, thermal units are recognized as being superior to cal- 
endar days for predicting dates of flowering and physiological maturity 
of corn in research and demonstration plots. However, in the realm of 
crop production forecasting at the regional or national level, one needs 
to know more than the rate of development of a specific corn genotype. 
He needs information about the status of the whole corn crop over large 
areas that may have many different planting dates, genotypes, and man- 
agement practices. The timeliness and reliability of this information 
influence many decisions of economic importance to individuals involved 
in producing, storing, marketing, or consuming corn products. 

The Statistical Reporting Service (SRS) of USDA acquires, summar- 
izes, and reports data on the progress of crops in each state at weekly 
or monthly intervals throughout the growing season. Additional informa- 
tion could be obtained by using daily meterological data and reliable 
models of crop development to assess the status of the crop in the 
region of interest. These models could be updated as needed using the 
data reported by SRS. However, the validity of using models of crop 
development for large areas has not been demonstrated. 

The objective of this research was to evaluate the use of statisti- 
cal data from SRS for assessing the development of corn in crop report- 
ing districts (CRD) of Indiana and Iowa. These data from SRS repre- 
sented means of adapted genotypes of corn in each CRD. Preliminary 
evaluation of the crop models used data acquired from research plots at 
an agricultural experiment station. 


MATERIALS AND METHODS 
Agronomy Farm 

Agronomic and meteorological data used in the first phase of this analy- 
sis were acquired at the Purdue University Agronomy Farm in 1979, 1980, 
and 1981. An adapted hybrid, Becks 65X, was grown on Chalmers silt loam 
(fine-loamy, mixed, mesic Typic Argiaquoll) at three densities (25,000, 
50,000, and 75,000 plants/ha) in 76-cm rows. Planting dates were 2, 1 6 , 
and 30 May 1979, 7, 16, 22, and 29 May and 11 June 1980, and 8 and 29 
May and 11 June 1981. Prior to planting, 200, 50, and 95 kg of N, P, 
and K per hectare, respectively, were applied. Stages of development 
(9) were observed once a week in 1979 and twice weekly in 1980 and 1981. 
Dates of silking and physiological maturity (black layer) were recorded 
when at least half of the plants of each planting date reached a parti- 
cular stage of development. 
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Daily meteorological data were recorded at the cooperative National 
Weather Service station (West Lafayette 6 NW) which was within 300 m of 
the plots. Daily maximum and minimum air temperatures were measured in 
a standard Cotton Region shelter. 


Crop Reporting Districts 

The percentages of the acreages planted, silked, and mature in each 
of the nine crop reporting districts (CRD) of Indiana were taken from 
the Annual Crop and Livestock Summary (16). Similar data for the nine 
CRD of Iowa were extracted from the annual Iowa Crops Weather Summary 
(17). Dates on which 25, 50, and 75? of the crop in each CRD reached 
each stage of development were linearly interpolated from these data 
(16,17). 

Meteorological data consisting of daily maximum and minimum air 
temperatures for 1969 to 1980 were selected for five National Weather 
Service (NWS) cooperative stations in each CRD of Indiana and Iowa (18). 
Stations with similar times of observation were selected to reduce any 
bias. Mean daily maximum and minimum temperatures for a CRD were com- 
puted from daily maximum and minimum temperatures reported by the five 
NWS stations in each CRD. The 12 years of data were assumed to repre- 
sent a random selection of years for each location and were divided into 
calibration (1969, 1971, 1973, 1975, and 1977) and test (1970, 1972, 
1974, 1976, 1978, 1979, and 1980) sets. 


Models and Analyses 

Four thermal indexing methods and the number of calendar days after 
planting (DAP) were evaluated for precision and accuracy. The first 
index, Growing Degree Unit (GDU), is the simplest thermal method and is 
defined as the daily mean air temperature minus a base temperature for 
growth of 10° C. The daily values of GDU are summed from the beginning 
to the end of each stage of development. For daily mean temperatures 
less than 10° C, GDU = 0. The dates that 25, 50, and 75? of the corn 
acreage had been planted in each CRD of Indiana and Iowa were used to 
start the accumulations of the thermal indexes. Dates that 25, 50, and 
75? of the corn acreage in each CRD had silked or reached physiological 
maturity were the ending dates. 

Modified Growing Degree Unit (MGDU) index (4) is the same equation 
as GDU but with a threshold of 30° C imposed on maximum temperature and 
a threshold of 10° C imposed on minimum temperature. 

Heat Stress (HS) index (7) is the same equation as MGDU but with a 
decrease in thermal unit accumulations for maximum temperatures greater 
than 30° C. 
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Function of Temperature (FT) index (5) is the mean of the relative 
growth rates for the daily maximum and minimum air temperatures. Four 
line segments which define FT are as follows: 


FT = 

0.027T - 0.162; 

if 

6° C 

< T 

< 21° C, 

FT = 

0.086T - 1.41; 

if 

21° C 

< T 

< 2S 3 C 

FT = 

1.0 ; 

if 

28 

< T 

< 32° C 

FT = 

-0.083T + 3.67; 

if 

32 

< T 

< 44° C 

FT = 

0 for 6° C > T > 

44° 

C. 




Daily FT was calculated as mean of the FT for the maximum temperature 


and the FT for the minimum temperature (5). The FT values used in this 
research were computed using air temperatures only rather than the com- 
bination of soil and air temperatures (5). 


The average thermal units and the number of calendar days accumu- 
lated from planting to silking, planting to physiological maturity, and 
silking to physiological maturity were calculated for the calibration 
years and used to predict dates of silking and physiological maturity 
for test years. Accuracy was measured as absolute errors in days, that 
is, the predicted date of stage minus the actual date of stage . Bias 
was measured as errors in days for predicted minus actual dates. Multi- 
ple range tests were used to separate significant differences in bias 
and accuracy among the models. 


RESULTS AND DISCUSSION 
Agronomy Farm 

The means, standard deviations, and coefficients of variation (CV) 
for the five models evaluated at the Agronomy Farm are shown in Table 1 . 
The GDU model had the smallest CV and the calendar days model had the 
largest CV for planting to silking. All the thermal models depicted 
silking better than calendar days for the wide range of planting dates 
used in the three years at the Agronomy Farm. 

The corn hybrid grown at the Agronomy Farm did not reach physiolo- 
gical maturity (i.e., black layer) before frost when planted after 10 
June in 1980 or 1981. Thus the statistics in Table 1 for physiological 
maturity are based on fewer observations than for silking. Differences 
in CV among the models were very small for planting to physiological 
maturity. However, for the silking to physiological maturity interval 
CV for the calendar days model was much smaller than CV for the thermal 
models. This observation is supported by Shaw and Thom (13) who noted 
that the interval from silking to physiological maturity is relatively 
constant over years. 
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Table 1. 
(CV) 

Means , standard deviations 
of thermal and calendar days 

(s), and 
models at 

coefficients of 
Purdue Agronomy 

variation 

Farm. 



Thermal 

. Models 


Calendar 

Statistic 

GDU 

MGDU 

HS 

FT 

Days 




Planting to Silking 



Mean 

818 

804 781 

37.2 

68.6 

s 

43 

45 50 

2.4 

6.4 

cv,? 

(n=1 1 ) 

5.3 

5.6 6.4 

6.4 

9.3 


Planting to Physiological Maturity 


Mean 

1499 

1497 

1466 

70.0 

133.0 

s 

70 

62 

60 

2.9 

6.0 

CV,? 

4.7 

4.2 

4.1 

4.1 

4.5 

(n=9) 








Silking to Physiological 

Maturity 


Mean 

676 

686 

677 

32.4 

62.7 

s 

64 

50 

44 

2.4 

2.6 

CV,? 

9.5 

7.3 

6.6 

7.4 

4.1 

(n=9) 


- 





Comparing models solely on the basis of CV of accumulated units for 
a number of environments provides an incomplete evaluation. A better 
way is to use the mean cumulative units from Table 1 for the respective 
models to predict the dates of silking and physiological maturity. Mean 
errors and mean absolute errors in number of days for the predicted date 
minus the actual date of each stage provide more realistic evaluations 
than simply CV. Mean error (e) is a measure of the bias of a model's 
predictions while mean absolute error del) measures its accuracy. The 
standard deviation of the absolute error (si~i) provides a measure of 
the precision or variability of a model's errors in predicting dates of 
corn silking or physiological maturity. Low variability signifies high 
precision. When silking dates of corn grown at Agronomy Farm were pred- 
icted, the thermal models were significantly more accurate than calendar 
days (Table 2). There were no significant differences among the thermal 
models. Rounding to the nearest whole day probably accounts for the 
slight positive bias (i.e., less than 1.0 day) exhibited by all of the 
models. 
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Table 2. Errors in days for predicted minus actual dates of silking at 
Purdue Agronomy Farm. 





Thermal Model 




Planting 




Calendar 

Year 

Date 

GDU 

MGDU HS 

FT 

Days 


Days 


1979 

2 May 

-5 

' -6 

-7 

-7 

-11 


16 May 

-1 

-2 

-3 

-6 

-2 


30 May 

-4 

-4 

-5 

-5 

3 

1980 

7 May 

2 

2 

3 

2 

-6 


16 May 

1 

2 

2 

2 

-1 


22 May 

1 

2 

3 

2 

2 


29 May 

0 

0 

1 

1 

2 


1 1 June 

-1 

0 

2 

1 

7. 

1981 

8 May 

5 

3 

3 

3 

-5 


29 May 

5 

5 

5 

5 

7 


11 June 

6 

6 

5 

6 

9 


e+ 

0.8a| 

0.7a 

0.8a 

0.4a 

0.5a 


lelt 

2.8b 

2.9b 

3.6b 

3.6b 

5.0a 


3 | - i t 

1 el 

2.2 

2.1 

1.8 

2.2 

3-3 


+ Mean error (e), mean absolute error ( I e I ) and standard deviation of 
mean absolute error (s|_|). 

| Within each line, means followed by the same letter are not signifi- 
cantly different at a = 0.05 level using Duncan's multiple range 
test. 


The errors and absolute errors for predicting physiological 
maturity dates using thermal models (Table 3) are at least double the 
errors for predicting silking dates using the same models (Table 2). 
The major source of the variation unaccounted for by thermal models in 
predicting date of physiological maturity appears to occur between silk- 
ing and physiological maturity. In contrast, the absolute errors for 
the calendar days model remain relatively constant for both stages of 
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Table 3. Errors in days for predicted minus actual dates of physiological 
maturity at Purdue Agronomy Farm. 


Year 

Planting 

Date 

Thermal 

Models 



Mixed 

Model s+ 


Cal. 

Days 

GDU 

MGDU 

HS 

FT 

GDU' 

MGDU' 

HS' 

FT' 




















- 



1979 

2 May 

-13 

-14 

-17 

-14 

-5 

-6 

-7 

-7 

-12 


1 6 May 

3 

0 

-3 

1 

3 

2 

1 

-2 

2 


30 May 

14 

11 

2 

7 

-4 

-4 

-5 

-5 

3 

1980 

7 May 

-6 

-4 

-2 

0 

5 

5 

6 

5 

-3 


1 6 May 

-6 

-2 

-1 

-2 

5 

6 

6 

6 

3 


22 May 

-6 

-2 

1 

-2 

4 - 

5 

6 

5 

5 


29 May 

-2 

4 

6 

4 

3 

3 

4 

4 

5 

1981 

8 May 

23 

10 

6 

9 

3 

1 

1 

1 

-7 


29 May 

26 

26 

24 

26 

2 

2 

2 

2 

4 


4 

3-7a § 

3.2a 

2.0b 

3.2a 

1.8b 

1.6b 

1.6b 

1.0c 

O.Od 


lel 

11.0a 

8.1b 

7. 1b 

7.2b 

3.8c 

3.8c 

4.2c 

4.1c 

4.9c 


3 |e| 

8.7 

8.2 

8.7 

8.4 

1.1 

1.9 

2.3 

2.0 

3.1 


+ The mixed models. predict physiological maturity by using thermal models 
to estimate date of silking and then adding the mean number of days from 
silking to physiological maturity from Table 1 . 

| Mean error (e), mean absolute error (|e|), and standard deviation of 

mean absolute error (s|-|). 

i e i 

§ Within each line, means followed by the same letter are not signifi- 
cantly different at a = 0.05 level using Duncan's multiple range test. 


development. The large positive errors observed in 1981 for all of the 
thermal models (Table 3) were due to very slow accumulations of thermal 
units late in the fall. Calendar days, on the other hand, accumulate 
uniformly. 
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A "mixed" model could exploit both the advantages of the thermal 
models for predicting silking dates and the reliability of the calendar 
days model for predicting physiological maturity. To test this mixed 
model concept, the original thermal models were used to predict silking 
dates and then the mean interval in days from silking to physiological 
maturity (from Table 1) was added to predict dates of physiological 
maturity. For example, the expected silking date would occur when 818 
GDU had accumulated after planting and the expected physiological matur- 
ity date would occur 63 days later (Table 1). 

A prime (') distinguishes the mixed models (thermal + days) from the 
conventional thermal models (Table 3). The accuracies of the mixed 
models are better than the accuracies of the conventional thermal or 
calendar days models for predicting physiological maturity of corn. 
There were no significant differences in accuracy among the mixed 
models. The mixed models appear to capitalize on the advantages of the 
thermal models for predicting date of silking and on the advantages of 
calendar days for predicting physiological maturity. 


Crop Reporting Districts 

The means of thermal and calendar days models for each CRD in Indi- 
ana are presented in Table 4 for the five calibration years. The number 
of thermal units accumulated for each interval increased from northern 
to southern CRDs while the number of calendar days remained nearly cons- 
tant or decreased slightly. Similar trends were observed in the data 
for Iowa and only the state means are presented in Table 5. 

Calendar days consistently had the lowest CV and GDU had the high- 
est CV for each of the three intervals in Indiana (Table 4) and Iowa 
(Table 5). during the calibration years. These results, using statisti- 
cal data from CRDs, contrasted sharply with our data from the Agronomy 
Farm (Table 1) and with many previous reports which have concluded that 
thermal units are significantly superior to calendar days in predicting 
dates of flowering (3,5,7,8,12,19). However, the trends for the inter- 
vals from planting to physiological maturity and from silking to phy- 
siological maturity (Tables 4 and 5) were consistent with trends 
observed at the Agronomy Farm (Table 1). 

One possible source of error introduced by using these statistical 
data for CRD was that the first 25 % of corn planted was assumed to be 
the first 25^ to reach all other stages of development. This assumption 
should be reasonable unless most farmers in a CRD shift to short-season 
corn genotypes as planting progresses. Such a shift is most likely to 
occur only in years when planting is delayed much later than normal. 
Other factors, not present in controlled experiments, also may affect 
analysis of statistical data on crops over large areas. For example, 
soil productivity and level of management may vary greatly from location 
to location and cannot be controlled by the investigator. The statisti- 
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Table 4. Means of thermal and calendar days models for planting to 
silking, planting to physiological maturity, and silking to physiol- 
ogical maturity of corn in crop reporting districts (CRD) of Indi- 
ana in calibration years. Data for CRD are means of three planting 
dates per year and 5 years (n=15). 


CRD 


Thermal 

Models 


Calendar 

Days 

GDU 

M3DU 

HS 

FT 



Planting to Silking 


NW 

728 

728 

688 

35.4 

67.6 

NC 

728 

72 8 

691 

35.5 

66.2 

NE 

717 

721 

693 

35.4 

66.3 

WC 

764 

759 

722 

36.8 

66.9 

C 

747 

747 

718 

36.5 

67.3 

EC 

735 

736 

708 

36.2 

67 .0 

SW 

840 

815 

775 

38.5 

64.8 

SC 

789 

777 

743 

37.6 

65.1 

SE 

850 

824 

785 

39.4 

65.7 

Mean 

767 

759 

725 

36.8 

66.3 

3 

73 

61 

62 

2.8 

4.0 

d,% 

9.5 

8.0 

8.6 

7.6 

6.0 



Planting to Physiological Maturity 


NW 

1385 

13 80 

1326 

67.6 

124.8 

NC 

1353 

1352 

1299 

66.6 

122.8 

NE 

1303 

1308 

1265 

64.9 

121 .2 

WC 

1452 

1440 

1386 

• 70.4 

125.1 

C 

1429 

1424 

1381 

70.1 

126.1 

EC 

1378 

1379 

1333 

68.1 

125.4 

SW 

1607 

1562 

1498 

74.3 

121 .9 

SC 

1461 

1435 

1382 

69.7 

118.7 

SE 

1563 

1516 

1451 

72.6 

120.4 

Mean 

1437 

14 22 

1369 

69.4 

122.9 

s 

113 

95 

91 

3.9 

6.1 

CV,% 

7.9 

6.7 

6.6 

5.6 

4.9 



Silking to 

Physiological Maturity 


NW 

657 

652 

638 

32.2 

57.2 

NC 

625 

624 

608 

31.1 

56.6 

NE 

585 

587 

572 

29.5 

54.9 

WC 

688 

681 

664 

33.6 

58.3 

C 

682 

677 

663 

33.6 

58.8 

EC 

643 

642 

626 

31.9 

58.4 

SW 

767 

747 

723 

35.8 

57.1 

SC 

671 

659 

639 

32.1 

53.5 

SE 

713 

692 

666 

33.2 

54.7 

Mean 

670 

662 

644 

32.6 

56.6 

s 

67 

60 

59 

2.7 

4.5 

OJ,% 

10.0 

9.1 

9.1 

8.2 

7.9 


+ Crop reporting districts are North West, North Central, North East, 
West Central, Central, East Central, South West, South Central and 
South East , respectively. 
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Tkble 5. Means, standard deviations (s) and coefficients of variation 
(CV) of thermal and calendar days models in Iowa during calibration 
years. Data are means of nine CRD, three planting dates, and 4 
years (n=108). 




Thermal Models 


Calendar 

Statistic 

GDU 

MGDU 

HS 

FT 

Days 



Planting to Silking 



Mean 

781 

773 

734 

37.6 

70. 1 

s 

59 

49 

50 

2.2 

3.0 

CV,* 

7.5 

6.3 

6.8 

5.8 

4.3 



Planting 

to Physiological Maturity 


Mean 

1376 

1355 

1293 

65.9 

120.8 

s 

79 

66 

67 

3.1 

5.2 

CV,* 

5.7 

4.9 

5.2 

4.8 

4.3 



Silking 

to Physiological Maturity 


Mean 

595 

583 

559 

28.3 

50.8 

s 

70 

58 

58 

2.5 

3-0 

CV,* 

9.0 

7.5 

7.9 

6.6 

4.2 

cal data 

on crops 

acquired by 

SRS represent the 

average 

genotypes , 

planting dates, soil 

productivity 

, and level of management of 

each CRD. 


When the models. were used to estimate dates of silking and physiol- 
ogical maturity for the calibration years in Indiana and Iowa, no signi- 
ficant differences were observed in bias or accuracy. However, evaluat- 
ing a model on the same data used to develop the model tests only the 
goodness of fit of the model to the original data and does not test the 
predictive ability of the model. For a more rigorous test, we assumed 
that years are random and divided the data into two series. The mean 
thermal units and calendar days accumulated during the calibration years 
for each CRD were used to predict the dates of silking and physiological 
maturity in 7 additional test years. 

The thermal models were significantly less biased and more accurate 
than the calendar days model for., predicting dates of silking in both 
Indiana and Iowa (Table 6) . Predicting silking date simply as the num- 
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Table 6. Mean errors (e), mean absolute errors (lit), and standard 
deviations of absolute errors (sigi) in days for predicted minus 
actual dates of silking in test years. Data are means of nine CRD, 
three planting dates per year, and 7 years for both Indiana and Iowa 
(n= 189). 


Location 

Statistic 

GDU 

Thermal 

MGDU 

Models 

HS 

FT 

Calendar 

Days 

Indiana 

e 

0. 5a+ 

0.0b 

— Days 

0.0b 

-0.8c 

-4. Id 


lei 

2.8b 

2.4c 

2.3c 

2.4c 

5.9a 


3,-1 

2.3 

2. 1 

2.0 

2.0 

4.3 

Iowa 

1 e 1 
e 

0.7a 

0.3b 

0.7a 

0.0b 

-2.5c 


lei 

3.6b 

3.5b 

3.3b 

3.4b 

4.7a 


3 |e| 

3.1 

2.8 

2.7 

2.5 

3.3 


+ Within each development stage and statistic, means followed by the 
same letter are not significantly different at the ct = 0.05 level 

using Duncan's Multiple range test. 


ber of days after planting produced a mean bias of -4.1 days in Indiana 
and -2.5 days in Iowa. The length of the average interval from planting 
to silking was slightly shorter in the test years compared to the cali- 
bration years. The biases of all the thermal models (except for FT 
model in Indiana) were positive and were within 0.8 days of the expected 
date. Thus mean air temperatures probably were slightly warmer for the 
planting-to-silking interval during the test years than during the cali- 
bration years. Differences among the thermal models were small. This 
contrasts with a previous report (5) which indicated that the FT model 
was clearly superior for predicting silking of corn. 

The calendar days model underestimated (i.e., negative bias) phy- 
siological maturity (Table 7) by approximately the same number of days 
as it underestimated silking (Table 6). The number of days from silking 
to physiological maturity changed little during calibration and test 
years. All of the thermal models had a positive bias for estimating 
physiological maturity (Table 7). The FT model was more accurate and 
less biased than other thermal or calendar days models. This extends 
the FT model -concept (5) to predict physiological maturity as well as 
silking,. 
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Table 7. Means errors (e), mean absolute errors (lei) and standard 
deviations of absolute errors (sig-i) of thermal, mixed, and calendar 
days models for predicted minus actual dates of physiological matur- 
ity in test years. Data are means of nine CRD, three planting dates 
per year, and 7 years (n=l89). 


Loc. 

Stat 


Thermal 

Models 



Mixed 

Models^ 


Cal. 

Days 

. GDU 

MGDU 

HS 

FT 

GDU' 

MGDU 

HS' 

FT' 























Ind . 

e 

7.4a{ 

4.0c 

5.4b 

2. 6d 

0. 8e 

0.3e 

0.3e 

-0. 5e 

-3.8f 


151 

11.0a 

7.3b 

8.2b 

5.9c 

3. 4d 

3.3d 

3. 2d 

3.3d 

6. 7bc 


S |e| 

8.6 

5.8 

7.4 

5.2 

2.7 

2.6 

2.6 

2.7 

4.5 

Iowa 

e 

4. 6a 

1.6b 

3.7a 

1.2b 

-0.4c 

-0.8c 

-0.4c 

-1.1c 

-3.7d 


lei 

11.0a 

8.0b 

9.8b 

7.8b 

5.3c 

5.2c 

4.9c 

4.9c 

6.4c 


S |i| 

9.4 

6.9 

9.1 

6.7 

4.0 

‘ 3.9 

3.6 

3.5 

4.7 


+ The mixed models predict maturity date in each CRD by using the ther- 
mal models to estimate silking date and then adding the mean number 
of days from silking to physiological maturity. 


| Within each line, means followed by the same letter are not signifi- 
cantly different at the a = 0.05 level using Duncan's multiple range 

test. 


Significant improvements in both bias and accuracy occurred when 
the mixed models were used to predict physiological maturity (Table 7). 
For example, the accuracy of the FT' model in Indiana was 3*3 days com- 
pared to 5.9 days for the conventional FT model. There were no signifi- 
cant differences among the mixed models. These results were consistent 
with our data from the Agronomy Farm (Table 3). The improved accuracies 
of the mixed models occurred mainly in years when the rate of accumula- 
tion of thermal units late in the season was much slower than normal. 
The mixed models predicted a date of physiological maturity whereas the 
thermal models accumulated the expected number of units for physiologi- 
cal maturity too slowly. 

This experiment evaluated thermal, calendar days, and mixed models 
to predict dates of silking and physiological maturity of corn. The 
results obtained using statistical data from CRDs were comparable to 
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those obtained using observations of plants in controlled experiments. 
In general, the data from CRDs may be used to extend and test models 
developed at agricultural experiment stations. 
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Table A1. f Dates (days from 1 January) of corn planting in Indiana. Table A2. Dates (days from 1 January) of corn silking in Indiana. 
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Table A3. Date, (day, from I January) of oorn maturity In Indiana. Table M. Date, (daya from 1 January) of corn planting In Iowa 
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Table A6. Dates (days from 1 January) of corn maturity in Iowa. Table A5- Dates (days from 1 January) of corn silking in Iowa. 
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